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Abstract 



^^ , The basic concepts of non-commutative probability theory are reviewed and appHed 

^^ ' to the large A^ limit of matrix models. We argue that this is the appropriate framework 

for constructing the master field in terms of which large N theories can be written. We 
Qj I explicitly construct the master field in a number of cases including QCD2. There we 



both give an explicit construction of the master gauge field and construct master loop 
operators as well. Most important we extend these techniques to deal with the general 



j^ ■ matrix model, in which the matrices do not have independent distributions and are 

coupled. We can thus construct the master field for any matrix model, in a well defined 
Hilbert space, generated by a collection of creation and annihilation operators — one 
for each matrix variable — satisfying the Cuntz algebra. We also discuss the equations 
of motion obeyed by the master field. 
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1 Introduction 

Theories invariant under U{N) (or 0(A^))|, in which the basic dynamical variables are A^^ 
dimensional matrices in the adjoint representation of the group, simplify greatly in the limit 
of large A^. In some cases the simplification is so great that the N = oo theory is solvable. 
Large A^ matrix models are of great interest for many reasons |]l|. First, QCD is such a 
theory, if we regard the number of colors as a free parameter. There is much evidence that 
the large A^ expansion of QCD correctly captures the essence of confinement and asymptotic 
freedom and that 1/3^ is a good expansion parameter. Second, matrix models have proved 
useful as devices for constructing string theories. Thus the control of the large A^ expansion 
of simple matrix models led some years ago to the non-perturbative solution of toy string 
models in dimensions less than or equal to two 0. In fact, QCD itself might be such an 
example, the large A^ expansion of QCD might be described by a string theory, a goal which 
has been realized in two dimensions [Q. Therefore it is important to explore and develop all 
available methods for controlling large A^ matrix models. 

One of the most appealing ideas to emerge in the study of the large A^ is that of the 
master field [§]. The idea is that there exists a particular classical matrix field such that 
the large A^ limit of all U{N) invariant Green's functions are given by their values at the 
master field. Thus the master field is analogous to the classical field, in terms of which all 
correlation functions are determined in the classical, h —>■ 0, limit; 1/A^^ playing the role 
of h. The argument for the existence of such a master field is simple. Consider a general 
matrix model. By this we mean a theory in which the dynamical variables are N x N 
dimensional Hermitian matrices Mj with an action S'[M(] that is invariant under the global 
U{N) transformation Mj — > UMiU"^ . Consider the correlation functions of U{N) invariant 
observables. If we have a denumerable set of variables then the most general invariant is a 
function of the traces of products of the matrices, i.e., O = -^Tr [Mi-^Mi^ . . . Mi^] (normalized 
so as to have a finite limit as A^ — > oo.) For a field theory with continuum fields, such as QCD, 



§In the large N limit there is no difference between SU{N){ov SO{N)) and U{N) (or 0{N)). 



we also consider continuous products of the matrix fields, such as Wilson loops. Denote the 
expectation value of O as 

(O) = Z^^ j\{dM,e-'^''^^^Ti [M,,M,, . . . M,J, (1.1) 

i 

where Z = JYlidMiexp{—S[Mi]). The important property of the large N limit is that the 
expectation value of a product of such invariant observables factorizes^, ||, §: 

(Oi02) = (Oi)(02)+0(l/iV2). (1.2) 

This can be proved in perturbation theory by the analysis of the Feynman graphs. In lattice 
QCD it can also be proved order by order in the strong coupling expansion. Consequently, 
the variance of any invariant observable vanishes in the large N limit, namely the probability 
that O differs from its expectation value is of order l/N"^ 

{{O - {Om = m - {{0)f = 0{llN'). (1.3) 

This must mean that the path integral measure is localized on a particular set of matrices — 
the master field — up to a U{N) transformation; just as in the classical limit the path integral 
measure is localized, infinitely sharply ash ^ 0, on the classical solution of the field equation. 
Given the master field, i.e., a set of "oo x cxd" matrices Mi, all the correlation functions 
of the invariant observables are then calculable as 

(0) = tr[M,,M,,...M,J, (1.4) 

where no functional integral need be done, we simply evaluate the trace of the product of 
master fields. 

In a gauge theory, in addition to the global U{N) symmetry that we used above we 
have a local U{N) gauge symmetry. In that case, when considering gauge invariant Green's 
functions, we can only conclude that the path integral is localized as A^ — i> oo on a single 
gauge orbit of the gauge group. In other words, if A^{x) is the master gauge field then an 
equivalent master field is A^{x) = U{x)A^{x)W{x) — iU{x)dfj,W{x). 



Given the master field for a pure gauge tfieory, say QCD in four dimensions, one could 
then calculate the meson spectrum very directly. Since in the large N limit the quarks play no 
dynamical role, quark loops being suppressed by 1/A^, the quarks are spectators and can be 
integrated out. Thus, for example, the meson propagator G{x, y) = {^ {x)T'^ (x)'^ {y)T'^ (y)) , 
where F is a matrix in flavor space, is given by 

Gix,y) = /l)(^^A,)e--^'^^^['^^^'^'™^""*"^"-H.,(a:)rv^(a;)^(,)rv^(,) 

= '^^^I^^^^JT^r^^-i— -b). (1.5) 

Thus, if we knew the master field, A^(x), we could calculate the meson spectrum for 
N = oo. In a gauge theory we can argue further that the master field can be chosen, by 
a choice of gauge, to be independent of space and time! This is reasonable if we think of 
the master field as the field configuration that yields the large A^ saddlepoint of the path 
integral. Since the action and measure are translationally invariant we might expect that 
the saddlepoint is translationally invariant, so that A^{x) and ^^(0) are equivalent up to a 
similarity transformation, 

A^{x) = e^^•%(0)e-^^■^ (1.6) 

where P^ plays the role of the momentum operator. If so, then we can perform a gauge 
transformation A^(a;) -^ A^{x) = U{x)A^{x)W (x) —iU{x)d^U'^ (x) , with U{x) = exp{iP-x), 
to derive a gauge equivalent, x^-independent master field 

4 = 4(0) + P^, F^, = [A^,A,]. (1.7) 

Thus the complete solution of large A^ QCD would be determined if we could write down 
four "oo X oo" matrices A^\ 

But what kind of matrix is A^? What does an "oo x oo" matrix mean? What is f/(oo)? 
To make these questions sharper let us consider a solvable example of a large A^ matrix model, 
a model of n independent Hermitian matrices, where the vacuum-to- vacuum amplitude is 
given by 

Z = j\{VM,e-''^'''^^''^\ (1.8) 



iv"(x) =/ ;^^^. (1,9) 



where V^j(Mj) is an arbitrary polynomial function of Mj. 

The case of n = 1 is the one-matrix-model, which is easily solved for N = oo . The 
invariant observables are class functions of M, determined by the eigenvalues rrii, i = 1 . . . N, 
which for A^ = oo yield a continuous function m{x = i/N). The matrix integral can be 
reduced to an integral over rrii, by diagonalizing M = QmQ^^ (where m = diag(?7ii, . . . , rrin)), 
and using the fact that T>M = PfiHi d'miA{mi)'^ , where T>Q is the invariant Haar measure 
on U{N) and A(mi) = Y\i<j{mi — rrij). The eigenvalue density, p{m) = ;^J^ is determined 
for N = oohy the saddlepoint equation, 

x-y 

For the simplest Gaussian potential, V{M) = |M^ the eigenvalues have the famous Wigner 

semi-circular distribution 

p{x) = — V4-x2. (1.10) 

zvr 

Thus, in the one matrix model we can say that the master matrix is an oo x oo matrix 
with eigenvalues rrii, where the rrii are determined by p(x), the solution of ( |1.9|) . If we now 
return to the n-matrix model, since the matrices are independent the eigenvalues of each are 
determined, so that we can say that the master matrices are 

Mi = Qim^^nl, where m'^^ = diag(mf\ m^*^ ...,m^§), (1.11) 

and the Qi are undetermined unitary matrices. These master matrices are perfectly adequate 
to calculate decoupled observables such as (tr Mf), in which the fij's do not appear. However 
the general invariant observable in this theory is a trace of an arbitrary product of different 
Mj's, namely (Or), where F denotes an arbitrary word, i.e., a free product of Mj's: 

Or = ^Tr[M,,M,,...M,,...] (1.12) 

Here the product does depend on the fij's and if we choose any particular fij's in ( p..llD 
we would not get the correct answer. Of course if we integrate over all values of the fij's 
with Haar measure then we get the right result, however this would not be a master field 
description. 



There is a direct, but rather ugly, way of deahng with this problem. Consider the case 
n = 2, with two independent matrices. Write each N x N dimensional master matrix as a 
block diagonal matrix oi K M x M-dimensional matrices, where N = M ■ K, 



Mi 
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1 . . . K are specific M x M unitary matrices chosen at random from the 



group (with Haar measure) and m^*^ = diag(m-[' 



« m« 
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rriM), with the 171^ determined 



(as M -^ oo) by the saddlepoint eigenvalue distribution. The expectation value of arbitrary 
words of Ml and M2 will now be correctly given by the trace of these master matrices when 
we take both K ^ 00 and M ^ 00. For example consider 



l(Tr[MiM2MiM2])= lim ^Tr;v[MiM2MiM2] 



1 ^ 

lim Y. TrM[V;m(^)l//m(2) V,.m(^V/m 



(2)1 



M,K^oo MK 
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j=i 



)(l)o(2)t 



where Vj = fl, 0,^ are a set of unitary matrices chosen at random from U{N). In the limit 



00 j^ Ef=iiv,)aAv^Uiv,),Avl)g,H, 



oi K ^ 00, the average of the product of the V's, lim;^_^o„ j^ ^'j=l\- j/a,uy j jc,a\ - jjc,j x- j 
is equal to the integral of the product over Haar measure, which for M ^ cxd is 



-t^ 



1 



T^V{Vj)a,b{Vj)c,d{Vj)ej{Vj)g^h = -T-j^[SadSbcSehSfg + SahSbg^ce^df] " JT^SadSbg^ehSfc- (1-15) 



This is simply the law of large numbers. Inserting this into ( |1 . 14| ) yields 
^Tr [MiMaMiMa])— (tt[M^] (ir [M2])+(Tr [Mi]) ^Tr [M|] - ^ (ir [Mi])'(Tr [M, 



N 



2 

1.16) 



which is the correct answer. Thus, the master matrices Mi,i = 1,2, given by ( |1.13| ), will 
yield all invariant Green's functions, i.e., arbitrary words made out of Mi and M2. This 
construction can be generalized to the case of an arbitrary number of independent matrices, 
at the price of imbedding the diagonal matrices m*^*^ in larger and larger block matrices. This 



is a very awkward construction. It indicates however the nature of the "oo x oo" matrices 
that will be required to represent the master field. 

Recently I. Singer [^ has presented an abstract existence proof for the master field for 
QCD2 and pointed out the relationship to the work of Voiculescu on non-commutative prob- 
ability theory]^. Indeed, Voiculescu's methods yield a much more satisfactory framework 
for representing the master field for independent matrix models [Q . More important we have 
been able to generalize these methods to deal with the most general matrix model, including 
QCD in any dimension, thus yielding an explicit representation of the master field for any 
and all matrix models. We do not mean that all matrix models are solvable, but rather 
that we can define a well defined Hilbert space and a well defined trace operation in which 
the master field of any matrix model can be explicitly constructed, if one possesses enough 
information about the solution of the theory. Although this construction can be viewed as 
repackaging it seems that the language that we shall review and develop is very appropriate 
to the N = 00 theory and might lead to new methods for constructing the master field, or 
equivalently for solving the N = 00 theory. 

In Section 2 we discuss the general framework of non-commutative probability theory 
developed by Voiculescu, define the notion of free random variables and the construction of 
an appropriate Hilbert space in which the master fields of models of independent matrices 
can be constructed. We explore this construction for the most general such independent 
matrix model. We note that the generating functional introduced by Voiculescu in his 
construction of the representation of a free random variable has the interpretation of the 
generating functional of planar connected Green's functions. We also show that the master 
field can be regarded as the solution of a certain master field equation of motion. 

In Section 3 we consider the explicit construction of the master field for some particular 
solvable gauge theories. We first find an alternative, manifestly Hermitian form of the matrix 
field for independent Hermitian matrix models. We reformulate the master equations of 
motion in a form that is more useful. We then discuss the simplest gauge theory, the one- 
plaquette model, which undergoes a large-N phase transition as a function of coupling. Here 



we will find two master fields, one for each region of coupling. 

In Section 4 we turn from theories of independent matrices to the general case of coupled 
matrices. Based on our interpretation of the generating functional introduced by Voiculescu 
in his construction of the representation of a free random variable we give a graphical proof 
of the construction of the master field for independent matrices. This argument can then be 
extended to deal with more general matrix models. We show that the master field for any 
number of coupled matrices can be formulated within the same Hilbert space as before and 
give its explicit construction. That is, 

If we can solve a matrix model then we can write an explicit expression for 
the master field as an operator in a well defined Hilbert space, whose structure 
only depends on the number of matrix variables. 

Section 5 is devoted to the construction of the master field for QCD2. Here we shall give 
an explicit construction of the master field and show that we can choose a gauge in which it 
is spacetime independent. 

In Section 6 we discuss an alternate description of QCD2 in terms of loops. We con- 
struct master loop operators based on the observation that simple loops corresponded to 
free random variables and that any loop could be decomposed into words built out of simple 
loops. The simple structure of QCD2 is then a consequence of the fact that these form a 
multiplicative free family. We use these master loop fields to recover the master gauge field. 

Finally, in the last section we shall discuss some of the many directions of research that 
are suggested by this construction. 

2 Non- Commutative Probability Theory 



Voiculescu has introduced the concept of free random variables for non-commutative proba- 
bility theory, which seems to be the appropriate mathematical framework for constructing the 
master field. We shall start by reviewing this framework, with no pretense at mathematical 



rigor. For more details we refer the reader to [F 



2.1 Free Random Variables 



For ordinary commuting random variables the notion of independence is simple, namely 
the probability measure of the random variables Xi factorizes, /i(xi, . . . ,a;„) = I\ifi{xi). 
Consequently the expectation value of products of functions of the x^'s factorize 



{fl{Xl)f2{x2) ■ ■ ■ fn{Xn)) = J K^l, • • • , Xn) fl{Xi) f2{x2) • • • fn{Xn) 

= n/M^.)/.(^i) = (/i(^i))(/2(a:2))...(/n(a:„)). (2.1) 

i 

For non-commuting random variables this definition is much too strong. There is a weaker 
definition, that of free random variables which is conceptually analogous to independence, 
though completely non-commutative. 

A non-commutative probability space is called free if 

The expectation value of products of functions of the non-commuting variables 
Mi vanish if the expectation value of all the individual functions vanish: 

r//.(MiJ) = for all k=l, . . . , n-1 
(/i(M,J/2(M,J.../„(M,J) = 0if 

[ and Ik t ^k+i k=l, . . . , n 

(2.2) 

Note that in the above product the neighboring functions must be of different random vari- 
ables. 

This is a much weaker condition than the previous definition where, because of factor- 
ization, the expectation value vanishes if any of the individual expectation values vanish. 
Nonetheless, it is a powerful restriction on the non-commutative probability space, that is 
sufficient to express all expectation values of products of different variables in terms of the 
individual expectation values. This is shown by considering the product 

([/i(M.J - (/i(M,J)] [/2(M,J - (/2(M,J)] . . . [UMJ - (/„(M,J)]) = 0. (2.3) 



Expanding this product one can express the expectation value of a product of n functions 
in terms of expectation values of ra — 1, ra — 2 . . . functions. Iterating this procedure one can 
express the expectation value in terms of individual expectation values. 

The expectation values of free random variables are obviously not symmetric under the 
interchange of different non-commuting variables. However, it is a remarkable fact that if the 
variables are free, i.e., (|2.2|) is satisfied, then the expectation value is cyclically symmetric. 
This can be proved using the same strategy we just employed, namely ( p.2| ) can be used to 
inductively show that if the expectation value of 2, 3, ... n variables is cyclic then it follows 
that the same is true for the expectation value of n + 1 variables. For details see |^ . 

The advantage of this definition is that independent matrix models in the limit of A^ = 
oo are free non-commuting random variables. To see this we denote 

tr[/i(M,J/2(M,J.../„(M,J] = hm l(Tr [/i(M,J/2(M,J . . . /„(M,J]), (2.4) 

N—i-oo iV 

where the expectation value is taken with the measure given by ( |1.8| ). Assume that the 
individual /j's have vanishing expectation value, tr [/^(Mj^)] = 0,k = 1, . . . ,n and consider 
the Feynman diagrams that contribute to the product in perturbation theory. A given 
matrix Mj must be contracted, when we use Wick's theorem, with the same Mj appearing 
in another, non- neighboring, term. Contracting two Mj's will split the trace into a product 
of lower order traces that, when A^ = oo , factorize. Thus one can prove the above claim 
inductively. Of course in this case tr is manifestly cyclic, as it must be for free random 
variables. 

We can use the fact that independent matrix models describe free random variables to 
disentangle the expectation values of arbitrary words. Thus, using the above method we see 
that 

tr [M1M2M1M2] = 2tr [Mi]tr [(M2)^Mi] + 2tr [Msjtr [{MiyM2] - (tr [Mi])hi [(Ms)^] 

-(tr [M2])'tr [(Mi)2] - 4tr [Mi]tr [(M2)]tr [M1M2] + 4(tr [Mi])\tT [M2])' 

-(tr [Mi])2(tr [M2])' = tr [Mi]tr [(M2)'] + tr [Msjtr [(Mi)'] - (tr [Mi]f{tT [M2])^ (2.5) 



which agrees with ( 1.16 ). Therefore we see that the notion of free random variables auto- 

9 



matically captures the content of Haar measure for independent matrix variables in the hmit 

oiN = 00. 

2.2 The Hilbert Space Representation of Free Random Variables 

Given a collection of free random variables, {Mi},i = l,...n, the correlation functions 
(M"^ . . . M^'' . . .) are linear functionals on the free algebra generated by the Mj's. There 
exists a very general mathematical construction that associates elements of a C* algebra 
(with a positive linear functional (p defined on it), with operators on a Hilbert space with 
a distinguished unit vector |fi).0 In the case of matrix models of Hermitian or unitary 
matrices there is a natural involution operation — the adjoint, so that we wish to consider 
cases in which the above free algebra is actually a C* algebra. States on this Hilbert space 
are generated by 

|M,^...M,J^M,^...M,Jfi), (2.6) 

where Mj's are the operators that represent the Mj's. The inner product on this Hilbert 
space is defined via the linear functional (p 

{A\B) = ^(A^B), (2.7) 

where 1^4) and \B) are states of the form given in (|2.6| ). In particular 

{n\m^ . . . mjn) = {n\M,^ . . . m,„) = 0(m,, . . . m,j. (2.8) 

In the case of matrix models where our linear functionals are expectation values with 
respect to the measure nj'PMiexp[— V^(Mi)], together with the trace, we recognize that the 
above apparatus is the appropriate framework for constructing the master matrix operators. 
We see from the GNS construction that the required Hilbert space is huge — a Fock-like space 
consisting of states labeled by arbitrary words in the Mj's. This is in agreement with our 
discussion of the master field above where we argued that the Hilbert space would have to 
be very large. 

^ This is the Gelfand-Nainiark-Segal(GNS) construction. See [|o) 
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For a one-matrix model-involving the matrix M the space is actually quite simple and 
can be described by states labelled by |fi), \M) = M|fi), \M^) = M^\n), . . . |M") = M"|fi). 
However, for a matrix model with n independent matrices Mi the Fock space of words is 
isomorphic to the an arbitrary ordered tensor product of one matrix Hilbert spaces. Note 
that the order is important since MiMaMsl^) ¥" MiMi,M2\^)- 

An ordinary Fock space of totally symmetric or anti-symmetric states is generated by 
commuting or anti-commuting creation operators acting on the vacuum. We might try to 
construct the above Hilbert space in an analogous fashion, by creation operators dj, for each 
Mj, acting on the vacuum \Vt). However, since the words are all distinguishable we would 
have to use creation operators with no relations, i.e., there would be no relation between 
a[aj and aja\. This is indeed the case. As shown in |^ the above Hilbert space is identical 
to the Fock space constructed by acting on a vacuum state with creation operators d|, one 
for each Mj, and that Mi can be represented in terms of d| and its adjoint dj. Specifically 
the Fock space is spanned by the states 

(alj"Mal)"^^---(4rHfi), (2.9) 

where 

a,\VL) = 0, did] = 6ij. (2.10) 

This is not an ordinary Fock space. There are no additional relations between different dj's 
or different dJ's, or even for djdj, except for the one that follows from completeness 

Y,dla, = l-Pn = l-\n){n\. (2.11) 

i 

In the case of the one-matrix model this implies that [d,d'''] = Pq. 

This algebra of the dj's and the dps is called the Cuntz algebra. It can also be regarded 
as a deformation of the ordinary algebra of creation and annihilation operators. Indeed it is 
the q = case of the q-deformed algebra 

djdj — qdjdi = 5ij, (2-12) 
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an algebra that interpolates between bosons (for g = 1) and fermions (g = —1). The 
above space can be regarded as the Fock space we would use to describe the states of 
distinguishable particles, i.e., those satisfying Boltzmann statistics!^ Working in such as 
space is very different from working in ordinary bosonic Fock spaces. In some sense it is 
much more difficult, since we must remember the order in which the state was constructed. 
Thus simple operators in ordinary Fock space can become quite complicated here. For 
example the number operator in the case n = 1 is given by 

N =:-%:■.= t{a^ra\ (2.13) 

1 - ata ^^ 

and obeys the usual commutation relations with a and with a^. The reason that even such 
a simple operator is of infinite order in a and d^ is that it must measure the presence of each 
particle in the state, thus it must be the sum of the operators {d'^)'^d^ that count whether 
a state has a particle in the k*"^ position. In the general case, for any n, the corresponding 
number operator is given by 

oo 

^ = E E dl---dlai^...di,. (2.14) 

fc=i n,...jfc 

Clearly we need to develop methods for working in such strange spaces. 

2.3 The Fock Space Representation of Mi 

It remains to show that we can construct an operator Mj, in terms of Sj and a} that reproduces 
the moments of the matrix Mj. Thus, suppressing the indices i, we wish to find an operator 
M{d, d^) in the Fock space so that 

tT[MP]= lim [vMe-^"^'^^^^^—Ti[MP] = (n\M{d,d^)\n). (2.15) 

Such an operator is clearly not unique, since we can always make a similarity transformation 
M — > S~^MS^ where S leaves the vacuum unchanged 5*1^) = |i7) and (i7|S'^"'^ = (i7|. 
II Greenberg has discussed such particles with "infinite statistics" [O 
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Voiculescu shows that we can always find such an operator in the form 

oo 

M{a,a)) = a + Y.^na^''^ (2-16) 

with an appropriate choice of the coefficients M„. To determine the coefficients we note that 



tr [M] = (fi|M|fi) = Mo; tr [M^] = (n\M'^\n) = Mi + M^, 
tr[MP] = (fi|MP|fi) = Mp + (polynomial in Mo, Ml,..., Mp_i). (2.17) 



Therefore we can recursively construct Mq, Mi, . . . , Mp in terms of tr [M], tr [M^], . . . , tr [M^]. 
To construct the explicit form of these coefficients we establish the following lemma. 
Lemma Given an operator of the form T = a + Y^'^q tnO^^ we associate the holomorphic 
function K = \ + EZoinZ"". Then 

n\F'{f)\n) = f^^^F[K{z)], (2.18) 

where C is a contour in the complex z plane around the origin. 

To prove the lemma it is sufficient to prove it for monomial F's, namely to prove that 
n (f]|f "'i|fi) = §c Y-K'^iz). But n (fi|f "-i|fi) = nTr [T^-^P^]. Then we use the fact that 



2m 

[T, a)] = [a, d"''] = Pn to write 



n (niT^-^in) = nTr [r"~^[r, a^] = Tr [T", a^], (2.19) 

where the last equality follows from the fact that Tr [T", d"''] = Y17=o r*[T, d"'^]T"~*~^ and the 
fact that [T", d^] is a trace class operator. Finally we use the fact that if Ty is the operator 
associated with the function f{z), that has the Laurent expansion / = J2'^=-oo fn^"" , i.e., 
ff = E:Li /-„a" + /o + E:Li /nS^", then 

TT[fj,f,] = j^^J{z)g\z). (2.20) 

It is sufficient to establish this formula for the case where f{z) and g{z) are monomials, then 
( p. 201 ) follows by additivity. Consider f{z) = /„2" so that Tf = /„d"''"'. Clearly Ti [Tf,Tg] 
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will vanish unless Tg = g_na"', i.e., g{z) = g-nZ "• Using d"'^"|m) = a)^'^^'^'""-'>\Vt) = \n-[- m), 

oo n— 1 



m=0 m=0 

-^ dz „dz- 



+ y^ [(m + n\m + n) — (m — n\m — n)] = —n = <p z"' — ; — . (2.21) 



Using this formula to evaluate ( |2.19| ) we establish ( |2.18| ) for polynomial functions, namely 



nin\T^-^n) = §ctrK''i^)- 



We now apply this formula to determine the form of the operator M that reproduces the 
moments of the matrix M. Assuming that we have found such an operator, so that ( p. 15 ) 



holds. Then we can express the resolvent, R{C)-, the generating functional of the moments 
R{0 ^ £ C— Hr [M-] = tr [-^] = J dx^, (2.22) 



n=0 

as 



RiO = - I ^H'iX) logic - A] = / ^^ = HiC). (2.24) 

Jc Zm Jc 2m C — X 



m = |C-"-' (n|M"|n) = |;^C-'-' i |^M"«(.) = - £ |i logic - M(.)l, 

(2.23) 
where M{z) = 1/z + J^MnZ"'. Now changing variables in the integral, M{z) = X, z = 
M-\X) =H{X) we have 

tc 2mC-X 

Therefore we find that M{z) is the inverse, with respect to composition, of the resolvent, i.e., 
R{M{z)) = M{R{z)) = z. 

This allows us to construct the master field for the one-matrix model explicitly, since the 
resolvent can be constructed algebraically in terms of the potential V{M). In the simplest 
case of a Gaussian, V{M) = 2^Tr [M^], we have 

G{z) = - — — = ^M{z) = - + az; M = a + aal (2.25) 

2a z + y/z^ — 4a z 

This form for the Gaussian master field can be made explicitly Hermitian by a simi- 
larity transformation, using the number operator constructed above. Indeed if we take 
S = exp[— ^ log aiV], then 

M -^ SMS-^ = v^[a + d^ = ^fax. (2.26) 
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2.4 Connected Green's Functions 

For a non-Gaussian one-matrix model the master matrix M = a + Z^^o ^n^^" will have an 
infinite number of non- vanishing M„'s. The function M{z) = l/z + J^MnZ" has, however, a 
simple interpretation. Let us recall the relation between the generating functional, G{j), of 
Green's functions and the generating functional of connected Green's functions , 

OO 1 oo oo 

G{j) = Ej"(tr[M"]) = -i?(l/j); Hj) = Ej"(tr[M"])conn. = Ej>"- (2-27) 

n=0 ^ n=0 n=0 

As shown by Brezin et.al. [^ the usual relation that tp = log[G] does not hold for planar 
graphs. Rather the full Green's functions can be obtained in terms of the connected ones by 
replacing the source j in ipiJ) by the solution of the implicit equation 

z{j)=Mz{j)). (2.28) 



Consequently, if one solves (|2.28|) for z(j) then 



GO) = ,H.O)) = ifl (1/.) ^ m/J) = 4,) - ^^^ = ^ = J- (2.29) 

J z{l/j) R{j) 

Therefore the the function iIj{z)/z is the inverse, with respect to convolution, of the resolvent 
R{z). But we established above that M{z) is the inverse of R{z). Consequently 

The master field function M{z) is such that zM(z) is the generating functional 
of connected Green's functions. 

This explains why in the Gaussian case zM{z) = 1 + az"^, since the only non-vanishing 
n-point function is the 2-point function, and why M{z) will be an infinite series in z for 
non-Gaussian distributions. Since the resolvent is a solution of an algebraic equation of 



finite order, for a polynomial potential, |12] it follows that M(z) is a solution of an algebraic 
equation as well. This interpretation suggests a direct graphical derivation of the form of 
the master field that we shall present in Section 6 and that will prove to be the basis for 
generalizing this construction to the case of dependent matrices. 
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2.5 Equations of Motion 

There are many ways in which independent matrix models can be solved. Saddle point equa- 
tions, orthogonal polynomials or Schwinger-Dyson equations of motion. The later approach 
is particularly simple and leads to equations of motion for our master fields. The Schwinger 
Dyson equations of motion for the one-matrix model follow form the identity 

d_ 
for an arbitrary function / (a sum of polynomials) of M. Using the fact that 



/^^E^]^{exp[-iVTry(M)]/(M),,} = 0, (2.30) 



Q n~l 

(M"),, = Y.iM%iM--^-%, (2.31) 



and the factorization theorem for iV = cxd , we derive for f{M) = M" 

(^Tr [V'iM)M-]) = X:(^Tr [M^])(^Tr [Pr'^~']). (2.32) 

These equations yield recursion relations for the moments of M that can be used to solve 
for the resolvent. 

The N = oo equations can be reformulated in terms of the master field as 

5 



{n\ 



V'{M) - 



dM 



f{M)\n) = 0, (2.33) 



for arbitrary f{M)- In this equation we must define what we mean by the derivative with 
respect to the master field. This is defined as 



6M ^--0 e 

so that 



{n\^^ ■ M"|fi) = J2{n\M^\n){n\M''~^-^\n). (2.35) 

6M j=o 



With this definition ( |2.33| ) is equivalent to ( |2.32| ). Below we shall recast these equations in 



a form that might prove more useful. 
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2.6 The Hopf equation 

The Hopf equation appears often in the treatment of large A^ matrix models. It arises in 
the collective field theory description of QCD2 |T9| , ^, where it determines the evolution of 
eigenvalue densities. It is also the equation of motion of the c = 1 matrix model [^ and 
governs the behavior of the Itzykson-Zuber integral [|I4[. We shall see that it arises very 



naturally in the context of non-commutative probability theory for families of free random 
variables. 

Let us first introduce the concept of an additive free family. Given two free random 
variables Mi and M2, with distributions /ii and /i2, their sum Mi + M2 has a distribu- 
tion /i3 denoted by /xi © /i2. A one parameter family of free random variables, such that 
l^ti ® /Ut2 = /^ti+t2 5 ^ill be called an additive free family. In ordinary probabilty theory the 
distribution of the sum of two random variables is given by the convolution of the two in- 
dividual distributions. However the Fourier transform is additive, i.e., we add the Fourier 
transforms of the individual distributions to get the fourier transform of the sum. The 
non-commutative analog of the Fourier transform isthe 7?.-transform that we have already 
encountered above. In section 2.3 we represented the free random variable M by the operator. 



M = a+Y, M„a^", (2.36) 

n=0 



with the associated series. 



1 °° 1 

Miz) = - + y MnZ^ = - + niz) (2.37) 

Then it is shown in p[ that 'R-{z) is additive Q Namely, if Mi and M2 are two free ran- 
dom variables with 7?.-transforms TZi and TZ2 respectively, then Mi + M2 has a distribution 
described by 

M = a + 7^(a"^), with 7^(^) = 7^l(2) + 7^2(^). (2.38) 

It immediately follows that for an additive free family, TZ{z) must be linear in t. Thus, 
**This also enables one to establish a central limit theorem for free random variables S. 
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for example, a free Gaussian additive family has 

M{z,t) = - + tz, (2.39) 

corresponding to the family of distributions / PMexp[— ^Tr M^]. In general, for an additive 
free family 

M{z,t) = - + tip{z). (2.40) 

where '^{^z) need not be linear in z. 

Consider the distribution for the free random variable A^(t) = N^^Mif) where A^o is free 
with respect to the M's which are Gaussian, but otherwise has some arbitrary distribution. 
Due to the additivity of 7^.(2;), 

7^7v(^) =7^o(-^)+t;^. (2.41) 

We shall show that the resolvent R{C,t), which is the inverse of N{z,t) = - + TZn{z) obeys 
the Hopf equation. 



To see this note that if 



then, 



C = - + 7^Jv(^) = - + 7^o(2)+t2 (2.43) 

z z 



1 flJ^ 

R{C,t) = R{- + no{z)+tz,t)=z => — |.= 
z at 

^ OR, dR , dC , OR OR dR BR ^^^^^ 

This explains the ubiquitous appearence of the Hopf equation in large A^ theories. In par- 
ticular we canunderstand the origin of the Hopf equation in the c = 1 matrix model ^. It 
is easy to see from this argument that if instead of being Gaussian M{t) were some other 
additive free family, as described by ( p.4(]| ), then the equation for the resolvent R{(, t) would 

be modified to 

dR , ,dR , , 

These are the collective field theory equations for these general families. 
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We will show in section 6.4 that the Hopf equation also arises in the case of multiplicative 
free families. This will explain why it appears in QCD2, where the Gaussian nature of the 
master field will be responsible for its occurence (though it will not be the resolvent that 
will obey the equation.) 

3 The One-Plaquette Model 

The master field representation that we have constructed for independent Hermitian ma- 
trices is not manifestly Hermitian. However, as we remarked, there are many equivalent 
representations of the master field. In this section we shall derive a manifestly Hermitian 
representation of the master field for independent Hermitian matrices and then apply this 
construction of the simplest model of unitary matrices, the one-plaquette model that exhibits 
a large-N phase transition |l^ . 



3.1 Hermitian Representation 

We shall now give a prescription, again not unique, to construct a Hermitian master matrix 
M{a,a^) = M^{a,a^) that reproduces the moments of the one-matrix model of Hermitian 
matrices. The idea is to express M as a function of the Hermitian operator x = a + a^ . But 
X represents the master field for a Gaussian matrix model. Thus writing M in terms of x 
is equivalent to expressing a matrix with an arbitrary distribution in terms of one with a 
Gaussian distribution. This can be done directly by a change of variables in the probability 
measure of M. 

Write the moments of the matrix distribution, given in terms of the density of eigenvalues, 
as 

tr [M"] = / rfAp(A)A" = [ ^v/4^^ A"(x), (3.46) 



where the function x(A) is a solution of the differential equation dx/dX = p{X)/y/A — x^. 
Therefore if we are given the eigenvalue distribution p(A) we can construct the master field 
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as 

d\(x) a/4 — x^ 

M = X{a + a^) = X{x) = M{x); where A(x) is determined by — ; — = — - — . (3.47) 

ax p{X) 

In the case of many independent matrices Mj, we can find the master fields in Hermitian 
form as Mj = Aj(xj) = Aj(aj + a|), with each Aj being determined separately from the 
distribution of eigenvalues of M,. 

The master fields in this representation also obey the master equations of motion dis- 
cussed above. It is amusing, and perhaps instructive for more complicated models, to refor- 
mulate these in a way that allows for the construction of the master field directly using the 
equations of motion. The equations of motion ( p.33| ) can be rewritten as 

{n\V\M{x))f{M{x)) - [n, /(M(x))]|fi) = 0, (3.48) 

where 11 will be defined to be the conjugate operator to M in the sense that 

[U,M] = Pn = \n){n\. (3.49) 

Note that on the right hand side of the commutator we have the vacuum projection operator 
and not the identity. Since M is Hermitian we can choose 11 to be anti-Hermitian . Thus 
in the case of the Gaussian potential, where M = x, we have 

M = a + d\ Il = p = -{d-d^). (3.50) 

With this definition we have that 

[n,/(M)] = -t-/(M). (3.51) 

oM 

Therefore the equations of motion are equivalent to 

[Q\V'iM{x))f{M{x)) - flf{M{x)) + f{M{x))fl\n] 
= (^n\V'{M{x))f{M{x))-2Ilf{M{x)) + f{M{x))Il\n) = 0. (3.52) 

But, since the states f{M{x))\Q) span the Fock space as we let / run over all functions of 
M{x), these equations are equivalent to the condition that 

V'{M{x)) - 2fl] \n) = 0. (3.53) 
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This equation can be use to solve for the master field, ie., given the potential V{M) solve 
( p.53| ) for a Hermitian operator M in the Fock space where 11 is conjugate to M. The first 
step, given an ansatz for M = M{x) is to derive an explicit representation of 11. To do this 
we first note that 

Xi — Xi 

where the labels on the x operators means that we are to expand the fraction in a power 
series in xi and Xr and order the operators so that all the xi's are to the left of all the Xr^s. 
Using this notation we can then write 

n = ^^-^iH^^p. (3.55) 

In this expression, when the operators are ordered, p appears to the right of all the x/'s and 
to the left of all the x^'s. 

To illustrate how this goes consider the Gaussian case where V'{M) = M. Take M = 
gix + g2X^ + ■ ■ ■ Then using ( p.55| ) IT = 1/ gip — g2/gi{xp + px) + • ■ ■ The equation of motion 



then reads 

[gix+g2x''-2/g^p+2g2/gl{xp+px)+- ■ ■]\n) = [{gi-l/gi)a^+{)g2-2g2/gl{d^f+g2+- ■ -W) = 0. 

(3.56) 
Consequently we deduce that gi = 1, g2 = 0, . . . ^ M = x, II = x. 

3.2 The One-Plaquette Model 

The one-plaquette model describes unitary matrices U with the distribution 

We shall derive a master field for U in the manifestly unitary form If = exp[iH{x)], where 
H{x) will be the master field for the eigenvalues of U, 



[^Tt [t/"]) = fdea{e)e'''^ = (fi|e^"^(*)|fi) (3.58) 
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The N = oo eigenvalue distribution was determined in [ITH] to be 



'■cos(^).M-sin2r^) A<2 



^(l + ^cos0) A>2 

Following the strategy described above we can construct H by the change of variables 



2-na{6)d6 = \^A — x'^dx and H{x) = 6{x). It immediately follows from ( p. 591) that for 
weak coupling the master unitary field is given by 



U = exp[2i sin"^ V Ax/8], for A < 2. (3.60) 

The phase transition is visible in the master field, since Xx/8 is a Gaussian variable, whose 
means square value exceeds one for A > 2, at which point U ceases to be unitary. In the 
strong coupling phase the master field is given by 

U = e'^^^\ where H{x) + ^ sin H{x) = -x V4 -x^ + 2 sin"^ - (3.61) 

2a 2 2 

This master field has the remarkable property that 

fi|[/"|fi) = 6n,0 + jiSn,l + <5n,-l). (3.62) 



4 The General Matrix Model 

So far we have discussed only independent matrix models where the action can be written as 
S = J2i Si{Mi) and there is no coupling between the various Mi's. We found that the master 
fields can be constructed in a Fock space in terms of creation di and annihilation operators 
a|, one for each degree of freedom, where the only relation satisfied by these operators is 
didj = 5ij. Now let us consider the most general matrix model with coupled matrices, for 
example QCD in four dimensions. One might think that it would be necessary to enlarge 
the Hilbert space in which the matrices are represented, or to modify its structure. This is 
not the case. We show below that we can construct the master field in the same space as 
before, with no new degrees of freedom or relations between the dj';s and a] 's. The only new 
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feature will be that Mj will be constructed out of all the Oj's and aj's, not just those with 
j = i. 

Let us go back to the construction of the master field for independent matrices and give 
a graphical proof that the master field defined by 



M. = a. + EC+'«^", 



(4.63) 



where zMi{z) = il)i{z) = 1 + J^'^Pi'z^ is the generating functional of connected Green's func- 
tions of the matrix Mi, i.e., ipii^z) = I]^o(^^ [^j"])'^") yields the correct Green's functions. 
Consider the most general Feynman graph that contributes to 



(lTr[M,,M,,M,3...Mj). 



(4.64) 



The most general contribution to such a Green's function can be drawn, as in Fig. 1, in terms 
of connected Green's functions. Fig.l represents a contribution to the N = oo Green's func- 
tion (Tr [Mj^M| Ml MIM4MIM4]), where the solid circles represent the connected Green's 
functions. We are using the standard double index line notation for the propagators of the 
matrices. 




Fig. 1 A contribution to (Tr [MfM^MiMlM4MlM4]). The sohd circles represent 
connected Green's functions. 
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What is special about these graphs is that none of the hnes cross, i.e. the points around 
the circle corresponding to the matrices Mj, in the order determined by the above word, are 
joined by lines that do not intersect. In that case the double index graph can be drawn on 
the plane and contains the maximum number of powers of A^. 

Now let us note that these graphs are in one-to-one correspondence with the terms in 
the expansion of (n\MfMlMiMlMiMlMi\n), with Mi given by ( ^^631) . Writing out the 



expression for this vacuum expectation value we find a contribution that exactly corresponds 
to the above graph, namely 

(n|(ai + ...)■ (ai + ...)■ (as + ...)■(... + ^^^ + ..•)■(••• + V'^af + ...)■ (as + •••) " 

•(... + i^lal + ...)■ (a4 + ...)■ (as + ...)■ (as + •••)■(••• + i^l^f + •••)• 

(. . . + i^la\ + ...)\n)= i^li^li^li^li^l (4.65) 



Conversely, every non-vanishing contribution to (fi|Mj^M^|MiM|M4M|M4|fi\ corresponds 
to a planar graph, such as that depicted in Fig.l. Start along the circle at a point denoted by 
(n|, then Ml must create a line labeled by 1. The next operator is Mi, so it can annihilate 
the line 1, with coefficient ipl, or create another 1 line as in Fig.l. Then the operator M2 
cannot annihilate the line 1, since 0,10,2 = 0. Then the next M2 must annihilate the line 2, 
with coefficient ip^-, as in Fig.l. Clearly there is no way we can get graphs with lines 1 and 
2 crossed, the Cuntz algebra ensures that this cannot happen. So on around the circle. 
This argument is a graphical proof that the master fields is given by ( [4.63 ), with ipi{z) 



being the generator of connected Green's functions of the Mj's. However, this argument does 
not really depend on the matrices being independent. Consider a general matrix model, with 
partition function 

Z = Jl[VM,e-^^'''^^'''-'^'\ (4.66) 

i 

where 5* will in general contain interactions between different matrices. The general graph 
that contributes to the Green's function has a similar decomposition in terms of connected 
Green's functions. An example of a contribution to (tr [M1M2M3M4M1M3M5M4]). is de- 
picted in Fig. 2 
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<ni 1 




Fig. 2 A contribution to (tr [M1M2M3M4M1M3M5M4]). The solid circles represent 

connected Green's functions. 
The only difference is that the solid circles, representing the connected Green's functions, 
now can involve matrices with different indices. We can construct a master field for each 
Mj, in terms of the same creation and annihilation operators as before, as long as we let Mj 
depend on all the creation operators. Thus if 



m = ai + j2^^. 



t,31,]2,-,jk^jl^J2 ■ ■ ■ ^jn^ 



(4.67) 



fc = l 



where 'ipi,ji,J2,-,jk i^ ^^^ connected Green's function: 



An,n,-,j, = ^(Tr [M,M,,M,, . . . M,J), 



(4.68) 



then the vacuum expectation values of products of Mj's will correctly reproduce the A^ = 
00 limit of the general matrix model. Note that the coefficients ipi,ji,j2,--;jk ^^^ cyclically 
symmetric in the indices. 

This argument shows that the master field exists as an operator in a Boltzmannian Fock 
space constructed with the use of a creation operator for each independent matrix field and 
gives an explicit expression for the master fields in terms of the solution of the theory. It 



25 



is not, of course, an explicit expression for the master fields-to do so would be to solve the 
theory. 

One new approach to solving large N theories might be to explore the equations of 
motion for the master field operators. In the case of coupled matrices the Schwinger-Dyson 
equations can also be formulated as 

'dS[Mi,M2,-.] S 



n\ 



DM, 5Mi 



/[Mi,M2,...]|fi =0, (4.69) 



where-r^ is defined as before and / is an arbitrary funtion of the Mj. We can also define 
a conjugate operator, IIj, to Mj as before, that satisfies [ni,Mj] = 6ijPn. If we can find a 
Hermitian representation of Mj then the above equations of motion are equivalent to 

'dS[M,,M2,...] 



-2n, 



\n) = 0. (4.70) 



It might be fruitful to explore this equation as a way of solving large A^ theories. 

5 Two Dimensional QCD 

We now turn to discuss the master field for the simplest gauge theory, two-dimensional QCD. 
There are two approaches to the master field that one might pursue. One is to construct a 
master loop field, Uc, that would be used directly to reproduce the expectation values of the 
Wilson loops for A^ = oo . We shall discuss this approach in the following section. The other 
approach is to construct directly the master gauge field Afj,{x). This is simple for QCD2, 
since in an appropriate gauge the theory is Gaussian and corresponds to an independent 
matrix model, albeit one with continuum labels. 

Consider two-dimensional Yang-Mills theory in fiat space. We can work either in Eu- 
clidean space, R2, or in Minkowski space M2. If we choose an axial gauge n^'A^ = 0, where 
n^ is a unit vector in any direction the theory becomes Gaussian. This is a legal gauge on 
i?2 or M2. We choose for convenience the gauge Ai{x) = and work in Euclidean space, in 
which case 

Fio{x) = F{x) = diAo; C = ^TiFixf = ^Tr {diAo{x)f . (5.71) 
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The field strength F is therefore an independent Gaussian matrix at each point in spacetime 
and we can easily construct a master field to represent it in the continuum Fock space defined 
by 

a{x)a){y) = 5'^{x — y)] a{x)\VL) = Q] \xi,X2., ■ ■ ■ iXn) = a){xi)a\x2) ■ ■ ■a\xn)\^). (5.72) 

In this space we can write 

F{x) = a{x) + a\x). (5.73) 

Alternatively we can work in momentum space, where a{p) = J ^ exp{ip ■ x)d{x) satisfies, 
a{p)a^{q) = 6'^{p — q). We can then solve the equation diAQ{x) = F{x) to construct the 
master gauge field 

Ao{x) = i-[a{x) + d\x)] = [ ^— \e-'P-^d{p) - e'^-^a^p)] . (5.74) 



di J 271 pi 

We can now make a gauge transformation to a gauge in which the master gauge field will 
be independent of x^. To this end we define the momentum operator, P^, in Fock space so 
that 

P^'\pi,P2,...,Pn) = [J2Pi] \Pl^P2,---,Pn)- (5.75) 

An explicit representation of P^ in terms of creation and annihilation operators is 

oo „ 

P^^ = P^^^ = Y^ d^p, . . . d^pkPta\vi) ■ ■ ■ a\pk)d{pk) ■ ■ ■ a(pi). (5.76) 

k=i-' 

Thus when the /c^'^term in the sum acts on an n particle state, for A; < n, it removes k 
particles from the state, measures the momentum of the k}'^ particle and then puts the k 
particles back in the original order. The momentum operator has the standard commutation 
relations with the creation and annihilation operators 

[p/^, d{p)] = -p^d{p), [P^, d\p)] = p^d\p) ^ [P^, Aq{x)] = -id^Aoix). (5.77) 

Therefore Aq{x) = exp{iP-x)AQ{x) exp{—iP-x) and, as discussed in the introduction, we can 
make a gauge transformation on the above master field, with gauge function U = exp{iP-x), 
to derive a spacetime independent master field: 

ii = A, Ao = Po + i(0) = Po + / ^- \aip) - d\p)] . (5.78) 
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In this gauge the field strength is given by 



2^ 



'a{p) + a\p)\ =Ao(0) 



(5.79) 



By means of a similarity transformation, we can rewrite Aq as 



Ao = Po + 



<Pp 
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a{p) + —a\p) 
pi 



(5.80) 



in which the second term we recognize as a sum of master fields for a continuum of Gaussian 
matrix variables, where the momentum space connected two-point function, the propagator, 
is 1/Pi. One has to be careful in using this field to introduce an infrared regulator for small 
Pi. This can be done by cutting out a small hole in momentum space or by a principle value 



prescription for the propagator. As explained in |T^, |I8| gauge invariant observables do not 
depend on the regularization. 

This master field satisfies the master equations [D^,F^v{x)] = S/SA,y{x). In the original 
axial gauge this means that 

6 



maiAoix) 



SAr,(x] 



■]/(^)|fi)=o. 



(5.81) 



where as before the operator derivative is defined as 



(5.82) 



SAo{x) 

These equations of motion can be used to show that the Wilson loop, which can be 
written in terms of the master field as, 

(Wc) = {n\T{exp[tg /' i%(t)]}|fi) = {Q\T{exp[tg /' A%x{t))xom\Q), (5.83) 

JO Jo 

satisfies the Migdal-Makeenko equations 0. Note that in the first integral, written in terms 
of the spacetime independent master field the path ordering is still necessary since, for non- 
straight paths, A^i^(t) do not commute for different t's. 



6 Master Loop Fields 
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6.1 Wilson loops 

One can alternatively describe the master field for QCD in terms of Wilson loops. These are 
manifestly gauge invariant and contain, in principle, all information about gauge invariant 
quantities. They are also the natural variables for a string theory formulation of QCD . 
One would therefore like to have master loop operators to describe the large N limit of these 
loops. 

From the point of view of the master field of QCD, the loop approach is, in general, quite 
unwieldy. The space of loops is too large and overcomplete. There is a lot of redundancy in 
defining master loop operators for every possible loop. The space of loops is much bigger 
than the space of points. It is hard to see what a 'basis' might be in this space. Moreover, to 
extract information about, say, the spectrum of meson bound states seems extremely difficult 
in practice starting from these loops (even in QCD2). 

Nevertheless, in QCD2, the loop space and loop variables have many nice and simplifying 
features (mainly due to the area preserving diffeomorphism symmetry of the theory). These 
features of the loop variables are not immediately apparent from the master connection that 
we constructed above. More importantly, they enable one to explicitly construct master loop 
operators that reproduce an arbitrary loop average fairly easily. Starting from these loops we 
can, by considering infinitesimal loops, recover the master field of Section 4. Alternatively, 
it will also be possible to start from the master field and derive the master loops without 
too much effort. 

6.2 Free Random Variables in the Loop Space of QCD2 

The main tool in trying to solve for Wilson loop averages are the Makeenko-Migdal loop 
equations 0. In 2 dimensions, they are especially tractable and Kazakov and Kostov have 
shown how the average of an arbitrary Wilson loop can be calculated for A^ = 00 using these 
equations. We shall approach the calculation of loop averages in a somewhat different and 
suggestive manner. We shall start by decomposing an arbitrary loop into a word built of 
simple loops, all originating at some common base point but with otherwise non-overlapping 
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interiors. (By a simple loop we shall henceforth mean a non-self intersecting loop on the 
plane.) We shall argue that these simple loops form a family of free random variables. This 
will enable us to calculate an arbitrary loop average in terms of < tr [f/^.] >'s where the 
Uc/s are the free random variables for simple loops Cj. Simple loops will thus form a basis 
in loop space, though they still contain too much information and are overcomplete. 

Let's first show that a set of simple loops, based at one point and non-overlapping, 
correspond to free random variables. For concreteness consider the loops Ci, C2, C3 based at 
the point P as in Fig.3. 

We denote the holonomies along the loops Ci by Ui, 



A»^ 



i.e ^- — < — -\J, < — ~^ 

Ui = P exp{i (f A^dx^") (6.84) ^^'p'^""^^ 

Note that the f/j's are U{N) matrices. We claim that p- a Three Simole Loods 

the Ui have independent distributions and hence are free random variables in the large A^ 
limit. One way to see this is to use the heat kernel action, which we know to be exact in the 
continuum limit, 

Z = JI[VUl UZp[Up]; Zp[Up]=Y.dnXR[Up]e-^'^''^^-, (6.85) 

L plaqucttcs R 

where the sum runs over all reresentations of U{N), xr is the character of the representation 
and C2{R) its second Casimir operator. We can choose a triangulation of the plane such 
that the given contours, Cj, are the borders of some of the triangles. The self similar nature 
of the heat kernel always allows us to choose such a triangulation. Then, when we come to 
use this measure to calculate averages of products of the Cj's, we can integrate out all the 
other link variables, (this is only true on the plane), leaving us with an equivalent measure 
^ = I Hi^^UiYli Z[Ui], where the product runs over all the simple loops, f/j. Therefore the 
resulting integrals over the Ui are over independent distributions. 

Naturally, this can also be seen directly from the loop equations. We shall give a rough 
sketch below. The loop equations are the N = 00 Schwinger-Dyson equations for Wilson 
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loops. 



d 



dx/, 5af,„ 



W{C) U=.ir)=l dx,{(j)S^^\x{(j) - x{t))W{C,,C2) (6.86) 



Here the L.H.S. refers to a variation of the area of the loop by Sa^i, at x = x(r). The 
R.H.S. vanishes unless x{t) is a point of self intersection in which case Ci and C2 are the 
two loops into which C breaks up at that point. (The / refers to the exclusion of a = r 
in the integral.) In the large A^ limit, W{Ci,C2) = W{Ci)W{C2) and, as shown by [^, 
because of the area preserving symmetry of QCD2, the loop equations simplify to (See Fig.4) 




+ 




-2 




-> = 



A^/ A 




Fig. 4 The Loop Equations for QCD2 



{dk + di-di- dj)W{C) = W{Ci)W{C2), {d, = 






(6.87) 



where the Ai are the areas that meet at the point P of self intersection, at which the loop 
splits up into Ci and C2. (the R.H.S. of Fig.4). The equations ( |6.871 ) form a closed set of 
equations that determine the loop average for a loop with n self intersections in terms of 
ones with a lesser number of self intersections. They can be solved recursively in terms of the 
loop average for a simple loop. The latter can be computed either from perturbation theory 
or by imposing appropriate boundary conditions on loops with a large number of turns. It 
has the value W{C) = e~^, where A is the area of the loop. 
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If we now consider two simple loops, as before, based at some point and with holonomies 
U and V, then it is possible to compute < tr [[/"i^/™! . . . lf"'kymkj > using the loop equations. 
They give a recursion relation for such a word of length 2k in terms of shorter words. We 
simply state the general expression. 

=< tr [U"'] >< tr [l/^^tr t/^^K""^ . . . U''^-'V"'^-'U''^V"''] > 

- < tr [f/"iV™i] >< tr [[/"2\/™2 ^ _ _ Jjnk-iymk-iunkymk^ ^ 
+ ...-< tr [f/"ll/™if/"2y™2 . . . Ijnk-iynik-i^ >^ ^^ pn.ym,^ ^ 

+ < tr [[/"ly^if/^av^^a . . . f/nfe-iymfe-if/nfc] ><- ^r [V""'] > . (6.88) 

This recursion relation can be obtained by combining the loop equations at the vertices 
where the loop on the L.H.S. breaks up into the loops represented by the various terms on 
the R.H.S. This relation is easily seen to be equivalent to a uniform (Haar) distribution (at 
large A^) for the relative angular integrals between U and V. Thus it implies that U and V 
are free random variables. Conversely, if we were to assume that the holonomies of simple 
loops are free random variables we could derive this relation from the defining property of 
such variables, as we discussed previously. 

The above recursion relation is a very useful expression that will enable us to calculate 
arbitrary loop averages rather efficiently. For k = 1, 

< tr [f/'^i\/"^i] >=< tr [f/'^i] >< tr [V^'] > (6.89) 

and for k = 2, 

KtTpVUV] >=< tr f/ >2< tr [V^2] > _ < trU >^< tiV >^ + < tr [U^] Xtr^ >2 

(6.90) 
which tallies with ( p.5|) . If we have more than two such simple loops then these relations 
can be applied repeatedly to reduce the average to a product of averages of powers of the 
individual [/,'s. 
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< tr [t/"] >=< tr [t/-"] >=< tr [t/^"] >= -L[„_iJnA)e-"^ = - / p,{l ' ^ 



The loop equations ( p.87| ) can also be used to compute < tr [[/"] > for a simple loop with 
holonomy U. Later we will obtain the same answer by other means as well. The answer can 
be expressed in terms of Laguerre polynomials L^ [|19], pO|l , 

1 <. ^~ 1 

(6.91) 
where A is the area of the loop and we have exhibited the integral representation for the 
Laguerre polynomials. The first few terms are displayed below. 

<trt/>=e-4, <tT[U^]>={l-A)e-^, < tip^] >= {1 - 3A + -A^)e-^ . (6.92) 

6.3 Decomposing a Loop Into a Word 

Having seen that simple, non overlapping loops, based at a point are free random variables, 
we now proceed to show how an arbitrary loop can be written as a word built out of such 
simple loops. In fact, for a loop with n self intersections, there are {n + 1) windows (i.e. 
enclosed interiors) and the word will be built out of {Ui}'s, z = 1, 2 . . . n + 1, which will be 
associated with these windows. In the interests of clarity and to avoid notational clutter, we 
shall illustrate the decomposition in a few representative cases and its general nature will 
then be apparent. 



The simplest self intersect- 
ing loop is the figure of eight 
in Fig. 5. The loops Ci and C2 
are simple loops and hence U ^ m- c ^ 




= 00 



and V are free random vari- 
ables. Therfore the loop average is simply 



Fig. 5 The Figure 8 Loop 



W{C) =<tT[UV] >=< tiU >< tiV >= e-^-^^^ . (6.93) 
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The first, non-trivial example is the loop in Fig. 6a. We shall introduce a notation for 
loops in terms of line i 

segments. Thus (13) 
denotes a and (31), b. 
The bar is to distin- 
guish it from (31) which 
will denote a~^ - the 
oppositely directed line 
segment to (13) and 
similarly (13) for b^^. 

Then, together with (12)=c and (21) = d, the loop itself can be written as 





Fig. 6a Overlapping Loops Fig. 6b Simple Loop Decomposition 



C= (13)(31)(12)(21) = (13)(31)(12)(21)(12)(21)(12)(21) 



(6.94) 



Where we have inserted (21)(12) = 1 and (12) (21) = 1 (two back tracking loops enclosing 
zero area). This is geometrically equivalent to the loop in Fig.6b. As a word we see that it 
is UV^ where U corresponds to the loop (13)(3T)(T2)(21) and V to the loop (12)(2T). These 
are simple loops and therefore, their holonomies U and V respectively, have independent 
distributions. Therefore, 



W{C)=<tT[UV'] >=<trf/><tr[\/'] >= e"-e-^^(l - A 



(6.95) 



which is the standard answer. 
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Fig. 7 A Loop and its Decomposition Into Simple Loops 

Next consider the contour in Fig. 7. With the points labelled as shown, 



C 



(12)(21)(13)(34)(45)(54)(41) 

(12) (2T) (14) (43) (31) (13) (34) (41) (13) (34) (45) (54) (41) 



Ui 



(12)(21)(14)(43)(31) (13) (34) (45) (54) (41) (14)(45)(54)(41) 



Ui 



U2 



Us 



(13) (34) (45) (54) (41) (14) (45) (54) (41) (14) (45) (54) (41) 



U2 

UIU2USU2U! 



Us 



U3 



(6.96) 



where we have again introduced backtracking loops so as to peel off successively, the loops 
corresponding to the different windows. Note that all these loops are based at the common 
point 1, because of the introduction of backtracking or 'thin' loops. Therefore the L^j's are 
free random variables and 



WiC) 



< tr [t/if/af/sf/af/l] >=< tr Ui >< tr [U2U3U2Ui] > 
= < tr t/i > (< tr U2 > V tr [f/|] > - < tr t/2 > V tr [f/3 >< tr U^] > 

+ < tr [U^] >< tr f/3 >< tr [f/|] >) (6.97) 

This can be compared with the loop average computed by usual means, once we express the 
moments of Ui in terms of the appropriate polynomials, using ( |6.91|) . 
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u, 




Fig. 8 A Non-Planar Graph Decomposed into Simple Loops 



Finally, consider a case which is 'non-planar' in the terminology of [|2T|, i.e., the loop 
depicted in Fig. 8. With the segments denoted as shown, 



C 



(12)(23)(31)(12)(23)(31) 

(12) (23) (32) (21) (12)(23)(31)(T2)(23)(32)(2T) (T2)(23)(3T) 



U2 



Ui 



(12)(23)(32)(21) f/2 (12)(21) (12)(23)(32)(21) (12)(21) f/i 



U3 

U3U2U7'U^'UiUi 



u: 



Us' 



Ui 



(6.98) 
(6.99) 



Once again, the f/j's are free random variables and therefore the loop average in this case is 

W{C) = < tr [U^U2U^^U^^UiUi] >=< tr [U3U2U^^U^^U4] >< tr f/i > 
= < tr f/i >< tr f/2 >< tr [U^^^U^'^UiUs] >=< tr Ui >< tr U2 > 

(< tr f/4"^ >< tr f/4 > - < tr U^;^ >< tr U^^ >< tr U4 >< tr f/3 > 



+ < tr f/3-^ >< tr f/3 >) = e-^^ie-^" - e-(^*+^«) + e"^^) (6.100) 



which once again reproduces the usual answer. 
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By now the general procedure must be apparent - we decompose the loop starting with 
the segments bordering the outside and form loops from each of the windows giving a co- 
herent orientation. It is posssible to characterise these words algorithmically in terms of the 
graphical structure, but that is not pertinent to our present purpose. We also see that this 
process of associating a word with a loop and then using the recursion relations (|6.88|) makes 
the computation of complicated loop averages rather simple, in fact, mechanical. 

6.4 The master loop operators 

We have decomposed an arbitrary loop with holonomy Ur into a word F built of simple,non- 
overlapping loops Ci and holonomies f/c. such that 

Ur = l[U^ (6.101) 

r 

Since the Ud are free random variables, we can associate master loop operators to them, t/c-. , 
by the general construction of Section 2. Then the master loop operator Ur that reproduces 
the loop average is 

W{C) = (q\ (ur = nf^ci) 1^) • (6-102) 

We shall now construct the loop operators U (supressing the contour labels) in the form 

oo 

f/ = a + ^ cjfc+ia^^ (6.103) 

fc=0 

The cjfc can be determined from the < tr [U"'] > which we saw, are given by ( |6.91| ). We claim 
that with 

a;, = (_i)^-i^A^-ie-^^/2 (6.104) 

we have 

(n|f/"|fi) = -Lf„_i)(nA)e-"^ =< tr [t/"] > . (6.105) 

We shall demonstrate this directly below. We can also represent If in an explicitly uni- 
tary manner, as in the one plaquette model. However, the manifestly unitary form is not 
particularly elegant and we shall not present it here. 
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Consider an infinitesimal rectangular loop (as in Fig. 9) 
of area AA = Ax At. The master loop operator associ- 
ated with it is ( |6.103|) for A -^ O.Tlius to lowest order 
in AA, we have 



U{AA) =a + (l 



AA. 



AAa\ 



(6.106) 




Fig. 9 An Infinitesimal Loop 

where a refers to the annihilation operator at the point x, d{x). We can equivalently represent 

this, by perfoming a similarity transformation, in the form 



U{AA) 



AA 



) + i\/AA{a + a) 



exp{iH); H = y/AA{a + a^ 



(6.107) 



This is the explicitly unitary form for U for the infinitesimal loop. Note that if we naively 
drop the term linear in AA, which arises from —^H"^ = —^AA{a + a^)"^ = ~^AA + ■ ■ ■, 



as being of higher order than the vAA, then U would not reproduce the correct leading 
behaviour in A A in < tr [f/"] >. 

But we also know that the holonomy U around such a loop is, in say, axial gauge 

U = {1 - iAo{x,t)At){l + iAo{x + Ax,t)At) = (l + idiAoAA) = exp{idiAoAA) . (6.108) 



Comparing (|6.108| ) and ( |6.107|) we have 



VAAdiAo = (a + d^) 



(6.109) 



This is equivalent to the master field of Section 4. Indeed, if we discretise the theory, i.e., 
smear the fields over plaquettes of size AA, then the action reads as 

1 



S 



Y: (AA)Tr (9iAo)l 



(6.110) 



plaquettes 



We see that yAAdiAo are Gaussian free random variables, represented by (d + d^), the 
result we obtained from the loop operator. Of course, this should come as no real surprise. 
It is somewhat less trivial to start from the master gauge field and to calculate the master 
loop operators explicitly for finite loops. In QCD2, we can do this rather easily since the 
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?7's have the special property of not just being free random variables but of also being a 
multiplicative free family. This is a concept analogous to the additive free family that we 
discussed in Sec. 2. 5. We shall briefly explain this concept. 

The product of two free random variables with distributions, /ii and ^2 is again a free 
random variable with some distribution /i3 denoted by /ii ® /i2- A one parameter family of 
free random variables, such that ^t^ ® fit2 = /^tit2? ^ill be called a multiplicative free family. 
(Or equivalently fig^ <S) yUs2 = /^si+s2) if we redefine the parameter t — > s = logt.) 

We claim that U{A) are a multiplicative free family with the area A playing the role 
of the parameter s. In other words, given two simple loops Ci and C2, based at a point 
and non-overlapping, f/c;^ ( A 1) [7^2(^2) has the same distribution as UcioC2{^i + ^2)- Again, 
there are many ways to see this. One is from the fact that the heat kernel action is self 
reproducing and exponentially dependent on the area of the plaquette. 

A more explicit way is to 
note that Uci{Ai)Uc2{A2) has 
the same distribution as does 
UcMiWhM^)- This is evi- 
dent from Fig. 10, where we 
see that since W has the same 
distribution as ^^^(^2) and both 
are independent of V. But VW'^ 
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— > 
. W ' 


' V 
— > 
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' V 

— > 


i 



Fig. 10 A Simple Relationship Between Loops 



is a simple loop of area equal to the sum of the two areas. 

This fact alone, actually enables us to construct the U solely from the knowledge of the 
master gauge field, i.e. from the knowledge of an infinitesimal loop. To do so we need a 
non-commutative analog of the Mellin Transform in ordinary probability theory, which is 
multiplicative for the product of two random variables. It turns out that one can define 



such a transform 



-the S-transform, such that S^^^Sf^^ 



S'^j,^^2- For a multiplicative free 



family Ssj^Ss2 = Ss-i^+s2 (dropping the /x's.) S{z) is therefore exponential in s in this case. 
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The function S{z) for a non-commutative random variable U is constructed as follows : If 

oo 

<p{j) = J2{n\u-\n)r, (6.111) 



ra=l 



then construct the inverse function x{^)> i-^- 4>{x{^)) = ^- The S-transform is defined as: 

S{z) = ^X{z). (6.112) 

Since the f/'s are multiplicative, we can use the S-transform of an infinitesimal loop 
to obtain the exact S-transform for one of finite area, knowing that it must necessarily 
exponentiate. For the infinitesimal rectangular loop in Fig. 9. we saw, in axial gauge, that 

U ={l + idiAoAA) = exp{idiAoAA). (6.113) 



Arguing backwards now, from the discretised action for QCD2, for which y/AAdiAo = 
(a -|- a''"), we have 

U = exp{iH) H = VAA{a + a^). (6.114) 

Now, 

00 00 

0(j) = Y: (fi|f/1fi) j" = E(1 - n'AA/2)r, (6.115) 



n=l n=l 

keeping only terms of order AA. The sum can be performed giving 



J__^ . (3-j ) 
1-J 2^(l-j) 



'/'(j) = -rh- —Jj^-TS^- (6-116) 



Equating this to z and solving for j = xi^) (again to lowest order in AA only) gives 

Xaa{z) = -^—{1 + AA{1 + 2z)). (6.117) 

1 + z 

Therefore SAAiz), defined as Xaa{z)^-^, is equal to (1 -|- AA{1 + 2z)) for an infinitesimal 
loop. For finite area, this exponentiates as expected, to give 

SAz) = e4(^+^^) ^ xa{z) = -^e4(^+2^). (6.118) 
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We can now use the Sa{z), which we obtained from the master gauge field, to give an 
alternative derivation of (|6.91| ) for < tr [f/"] >. Since 



{e~'') = ^ < tr [[/"] > e"'"' 

n=l 
1 ' ' ' -ie\Ane, 



=^<tr[f/"]> = — / 0(e-*')e*"^rf^, (6.119) 

we have, with x{.z) = e~*^, 

< tr [U^] >= j ^^z[x{z)]-i-+^)^{z)^ (6.120) 

or on integrating by parts 

<tr[f/^>=i/|^[x(^)]-^ (6.121) 

In the case at hand Xa{z) = ^-|-e 2'*^^+^^) and hence 

< tr [f/(A)"l >=-{ —(1 + l)"e-"4(i+2-) (6.122) 

n J 2ni z 



which is ( |6.91 ). 



We can now see why the Hopf equation arises in the QCD2. In fact it, or its generalization 
( p.45| ), will appear for any multiplicatively free family of random variables. Define the 



resolvent of U{A) to be 



i?((, A) = ^ < tr [f/(A)"] > (-("+!). (6.123) 

n=0 



By definition (/) is related to the resolvent as 

i?(C, A) = i(0(i) + 1) ^ R{^ A) = x{z){z + 1). (6.124) 

C C x{z) 

Now exp[i^] = l/x{z) = (1 + z)/zexp[—A/2{l + 2z)], from which it follows that 

e{z)=iA{z + ^)-tlog^^^. (6.125) 



Therefore, using (|6.124|) and the fact that (j){x(z)) = z, we have 

R{e'' = ^,A) = e-'\<P{e~'' = x{z)) + 1) = e-'\z + 1) ^ e''R{e'\ A) = l + z. (6.126) 
X\Z) 
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Redefining w = i{z + ^) we have 

e{w) = Aw -Hog ^^^^^. (6.127) 

Then F{9, A) = i[e'^R{e'^, A) - I] = w satisfies 

dF dF , , 



w 



We can also employ Sa{z) to explicitly compute U in the form given by (|6.103|) . Thus if 



-| oo 

U{z) = - + J2iOk+iz'' (6.129) 

^ fc=0 

then U{z) is the inverse of the resolvent and 

X{z) X{z) 

Therefore 

;7(^e4(i+2-)) = l±ie-#(i+2^). (6.131) 

z 

Consequently, if z{y) is determined from 

^e4(i+2-)=y, (6.132) 

then 

Uiy) = -il + z{y)). (6.133) 

This is essentially the result we need. We have obtained U{y), albeit thus far in an implicit 
form. To obtain the coefficients u^, we must examine the relation (|6.132 ) more carefully. 
We have 



yU{y) = l+ziy) = l + Y^u;ky' 

k=l 

oo oo 

= 1 + X] a;,z'^e'=#(i+2.) ^ ^ = ^ Ukz'e'^'-'+^'l (6.134) 

fc=l k=l 

This is what determines the coefficients Uk- In fact, if we redefine ujk = A^~^e~^~Ck and 
z -^ Az, then ( |6.134| ) becomes an equation for the c^'s. 

oo 

z = J2ckZ^e''\ (6.135) 

fc=i 
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All the area dependence is gone and the c^ are just some numbers determined recursively by 
the above equation. In fact, the recursion relation is non trivial 

Cfc = - ^ Ck-r j (6.136) 

It can be checked that the Ck are precisely (—1)^^^^-^ due to the non-trivial combinatorial 
relation 

(6.137) 






k\ [k — r)\r\ 

But we can actually, rather simply argue that u^ have to take this form. That is, it is simply 
the highest power of A term that appears in < iiU^ >. This follows since < trf/^ > is, 
a polynomial with highest power of A being A^~^ multiplying, of course, the e~^"^/^. Since 
uJn oc A"-^^, the only term in < 0|f/^|0 > that can contribute an A^"^ term is uJk (other 
polynomial terms are of the form uji-^uJi2 . . .Wi^ with J2r''-r = k and thus are always of lower 
order). Therefore Uk must be precisely the term appearing with A''~'^ in < trf/^' > which is 
exactly what was given earlier (as can be checked from the expansion of the LJ _^n). 

7 Conclusions 

In this paper we have reviewed the basic concepts of non- commutative probability theory 
and applied them to the large N limit of matrix models. We discussed at length the work 
of Voiculescu on the properties and representation of free random variables. Since indepen- 
dent matrix models at A^ = oo are free random variables this appears to be the appropriate 
framework for constructing the master field. We discussed some of these models, including 
the one-plaquette model where we explicitly constructed the master field. We also discussed, 
at length, QCD2. In an axial gauge this theory can be regarded as a theory of independent 
matrices and thus we could give an explicit construction of the master gauge field. We also 
showed that there exists a gauge in which the master gauge field is spacetime independent. 
We also constructed master loop operators based on the observation that simple loops cor- 
responded to free random variables and that any loop could be decomposed into words built 

43 



out of simple loops. The simple structure of QCD2 is then a consequence of fact that these 
form a multiplicative free family. 

The most suprising and exciting of our results, however, is the extension of these tech- 
niques to deal with the general matrix model, in which the matrices do not have independent 
distributions and are coupled. Based on our observation that the generating function, intro- 
duced by Voiculescu to construct the representation of an independent random variable, can 
be identified as the generating function of connected planar Green's functions, we were able 
to construct the master field for any and all matrix models. Remarkably the Hilbert space 
in which the master fields are represented is unchanged — it is the Fock space generated by a 
collection of creation and annihilation operators satisfying the Cuntz algebra — one for each 
matrix variable. 

From some points of view our construction is somewhat dissapointing. First, although 
we have an explicit construction of the master field for any matrix model in a well defined 
Hilbert space, to actually write the master field explicitly would require a knowledge of all 
the connected Green's functions, which is tantamount to solving the theory. Thus from 
this point of view all we have done is to repackage the unknown solution. Second, we have 
almost as many degrees of freedom as before. The Hilbert space in which the master field is 
represented is almost as big as the full Hilbert space of the quantum field theory-i.e., there 
is an independent creation operator for each independent field variable. The only reduction 
is by a factor of A^^, since the large A^ limit has been taken. 

However, we believe that this reformulation is valuable. Clearly this is the appropriate 
framework for formulating the N = 00 theory. It also suggests new approaches towards 
solving the theory by constructing the master field — now a well defined operator in a well 
defined space. For example, one approach might be to explore the operator equations of 
motion for the master field, as we have discussed above. Do these, for example, follow from 
some kind of variational principle that could be the basis for an approximation scheme? Can 
one develop similar techniques for the Hamiltonian formulation of large A^ theories? 
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